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Abstract. 

We initiate the study of Majorana fermion codes. These codes can be viewed as 
extensions of Kitaev's ID model of unpaired Majorana fermions in quantum wires 
to higher spatial dimensions and interacting fermions. The purpose of Majorana 
fermion codes (MFCs) is to protect quantum information against low-weight fermionic 
errors, that is, operators acting on sufficiently small subsets of fermionic modes. We 
examine to what extent MFCs can surpass qubit stabilizer codes in terms of their 
stability properties. A general construction of 2D MFCs is proposed which combines 
topological protection based on a macroscopic code distance with protection based on 
fermionic parity conservation. Finally, we use MFCs to show how to transform any 
qubit stabilizer code to a weakly self-dual CSS code. 
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1. Introduction 

The physical realization of systems with topological quantum order and their theoretical 
description have been a topic of significant attention lately. This attention is partly 
motivated by the potential use of topologically ordered systems as fault-tolerant 
hardware in a quantum computer [1] . Encoding of quantum information into the ground 
states of such systems permits exponential suppression of dephasing, while the presence 
of a gap above the ground state suppresses thermal error excitations. The original 
insights [H [2] concerning the zero-temperature stability of topologically ordered systems 
have been fully rigorously proved for quite general quantum spin systems in Ref . [3j H] . 

A more ambitious goal for robust quantum information processing is to genuinely 
store and manipulate quantum information at nonzero temperature T > 0. Whether 
this is physically possible is the topic of the discussion on 'self-correcting' quantum 
memories [21 EH [7], which do not need to be continuously error-corrected as the standard 
theory of quantum fault-tolerance prescribes [8]. The question of self-correction goes 
under alternative guises such as the question of thermal stability [9j [10] or thermal 
fragility [XT] of quantum memories or the persistence of topological order at finite- 
temperature [T2l[T3l[T3]. 

In understanding aspects of topological order or its possible extension to finite 
temperature, it is important to study physical 'toy' models such as the well-known 
surface code family or general quantum error-correcting codes with geometrically-local 
generators. These toy models both teach us what may be possible at the conceptual level 
as well as pose an interesting challenge to engineer these Hamiltonians at the physical 
level. 

In this paper we introduce a class of toy models that can be viewed as extensions of 
Kitaev's ID model of unpaired Majorana fermions in quantum wires [15] to higher spatial 
dimensions and to interacting fermions. These toy models which we call Majorana 
fermion codes can be described by fermionic term-wise commuting Hamiltonians 
composed of geometrically-local interactions on a D-dimensional lattice. The purpose of 
Majorana fermion codes is to protect quantum information against low-weight fermionic 
errors, i.e., operators acting only on sufficiently small subsets of fermionic modes. One 
distinction between fermionic systems and the systems composed of qubits or spins is 
the presence of superselection rules. In particular, if a fermionic system interacts with a 
bosonic environment, conservation of the parity of the total number of fermions restricts 
the set of physically realizable errors to the so-called even fermionic operators. 

One question addressed in the present paper is whether the superselection rules 
permit more robust storage of quantum information based on Majorana fermion codes 
as compared with qubit stabilizer codes under the same geometric locality constraints. 
We partially answer this question in the negative by generalizing the no-go theorem for 
quantum self-correction based on 2D stabilizer codes [7J to Majorana fermion codes. 

On the positive side, we construct interesting 2D generalizations of Kitaev's model 
of unpaired Majorana fermions in quantum wires [15]. Specifically, we construct a 
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family of Majorana fermion codes encoding one qubit into a 2D lattice of fermions 
folded into a cylinder. The corresponding Hamiltonian has exactly two zero-energy 
Majorana modes, i.e., odd fermionic operators Cq,C\ supported on the opposite edges 
of the lattice, commuting with the Hamiltonian, and anti-commuting with each other. 
In coding theory language, Co and C\ are the logical operators of the code. The main 
advantage of the 2D model is that the logical operators C , C\ have a macroscopic weight 
proportional to the radius of the cylinder. It endows the encoded qubit with an extra 
degree of protection related to the macroscopic distance of the code and which does not 
rely on superselection rules. By varying the radius and the length of the cylinder one 
can combine the two types of protection in a controllable way. 

The additional protection by the code distance is completely analogous to the 
protection of quantum information encoded into the ground state of topologically 
ordered systems as discussed in [H EJ [3j Hj. Such protection might be necessary for 
example if one tries to build a multi-qubit register composed of ID quantum wires with 
unpaired Majorana modes. In this case a perturbation can couple the unpaired modes 
that belong to adjacent wires without violating the superselection rules and hence an 
additional protection using a Majorana fermion code with a large distance could be 
helpful. In addition, the superselection rule prohibiting odd error operators is not likely 
to be completely rigorous, for instance, if the environment supports gapless fermionic 
modes that can couple to the system, or, when a single unpaired fermion could jump 
from the superconductor onto the topological insulator (although such processes are 
energetically suppressed by a gap at low-enough temperature). 

Finally, we argue that Majorana fermion codes can be used as a tool to prove new 
facts about qubit stabilizer codes. In particular, we prove that any stabilizer code can 
be locally mapped onto a weakly self-dual Calderbank-Shor-Steane (CSS) code. The 
mapping preserves all parameters of the code such as the number of physical and logical 
qubits, the distance, and locality of the generators up to a constant factor. 

An open problem which is not addressed in our paper is by what physical means 
and mechanisms (interacting) Majorana fermion codes can be realized. One expects 
that, similar as for spin-systems, such models could emerge as effective many-body 
Hamiltonians for interacting fermion systems which are treated with a perturbative 
or renormalization-group flow analysis. There is ample physical evidence that 2D 
strongly-correlated electron systems support topological order; the question is whether 
the Majorana fermion code framework can help in understanding how such topological 
order emerges from the basic physical interactions. 

The paper is organized as follows. Section [2] reviews Kitaev's ID model of unpaired 
Majorana fermions and provides some motivation behind the present work. Section [3] 
introduces notations and necessary facts from the theory of stabilizer codes. A formal 
definition of Majorana fermion codes is given in Section HI The mappings between 
stabilizer codes and Majorana fermion codes are described in Section [5j This section 
also proves the equivalence between general stabilizer codes and weakly self-dual CSS 
codes. In Section [6] we discuss quantum error correction in the presence of superselection 
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rules. The generalization of Kitaev's ID model with unpaired Majorana fermions is 
described in Section [7J We give more examples of Majorana fermion codes that possess 
odd logical operators in Section El 

2. Why Majorana fermion codes? 

In Ref. [15] Kitaev considered the toy Hamiltonian of a ID chain of spinless fermions 
interacting with a superconductor. The interaction with the superconductor allows for 
the creation and annihilation of pairs of fermions so that the total Hamiltonian of the 
system preserves only the parity of the number of fermions. Any Hamiltonian involving 
spin or spinless fermions can be written in terms of Majorana fermion operators by 
taking = \{c 2 h-\ + ic2k) and a\ = \{c 2 k-i — ic 2 k). The Hermitian Majorana operators 
obey the relations 

For a particular choice of couplings, Kitaev's Hamiltonian on L spinless fermions, hence 
2L Majorana fermions, reads 



i.e. all Majorana modes are 'paired', except the first mode C\ and the last c 2 l- Thus 
this quadratic fermion Hamiltonian commutes with the operators c\ and c 2 l which we 
identify as the logical operators of a protected qubit. Alternatively, it can be said that 
the symmetry gives rise to the presence of a pair of zero-energy Majorana boundary 
modes which lead to degeneracy at the Fermi-level. Recently, there have been various 
proposals to realize such ID model with zero-energy Majorana fermions, for example 
at the boundary between a topological insulator and a superconductor [16J or in a 
semiconducting heterostructure [T7] . 

Ground states of H are —1 eigenvectors for every term ic 2 jC 2 j + \ in Eq. (CQ). The 
two-fold degenerate ground subspace has a basis |0) and |1) which satisfies 



The logical Pauli operators for a qubit encoded into |0) and |1) can be chosen as X = c%, 
Y = —c 2 l, and Z = ic\c 2 L- Note that two of the logical operators C\ and c 2 l are 
of odd weight and hence would require the coherent creation/annihilation of a single 
fermion which is prohibited by superselection. An essential part of the model is that 
the only even-weight logical operator C\C 2 l is very non-local. It is natural to assume 
that elementary perturbations to the Hamiltonian and errors can be represented by 
local even weight Majorana fermion operators. Hence in a perturbative analysis such 
as the Schrieffer-Wolf perturbation theory, the first contributions that split the energy 
degeneracy between |0) and |1) are expected to occur in 0(L)th order implying that 
the splitting in degeneracy between |0) between |1) is exponentially small in L. 



L-l 




(1) 



^CiC 2 l|0) = |0) and %c x c 2L \\) 
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The spectrum and properties of Kitaev's model -as any other quadratic fermion 
Hamiltonian in the theory of topological insulators- are efficiently computable and 
quantum circuits which employ only non-interacting fermion Hamiltonians and simple 
fermionic measurements are efficiently simulatable classically, see [18J. Hence, if we 
are serious about using fermionic systems to robustly store and manipulate quantum 
information (see e.g. [IS]), we will need some source of interaction to obtain quantum 
universality (see e.g. [20J ) . Let us mention that generalizations of the Hamiltonian 
Eq. (p]) to interacting fermions have been recently considered by Fidkowski and 
Kitaev [21] to study the effect of interactions on the classification of ID topological 
insulators. 

The toy model, Eq. (JT]), demonstrates that fermionic parity conservation provides 
an alternative protection mechanism for the encoded qubit unrelated to topological 
quantum order. It is therefore natural to ask whether topological protection can be 
combined with protection based on superselection rules in an advantageous way. The 
Majorana fermion codes introduced in the present paper provide a natural framework 
in which such a hybrid protection can be studied. 

3. Definitions and notations 

We review a few standard definitions and notations which are used in this paper. Let 
Xi,Yi,Zi represent the three Pauli matrices on qubit i. Let V n = (il, X\, Z%, . . . , X n , Z n ) 
be the Pauli group on n qubits generated by single-qubit Pauli operators and the phase 
factors ±1, ±i. The support of a Pauli operator P, Supp(P) is the set of qubits on 
which it acts non-trivially. The size of the support, |Supp(P)|, is also sometimes called 
the weight of P, denoted as \P\. 

A pair of logical operators for an encoded qubit is denoted as (X, Z) and Y = iXZ. 
A stabilizer code is determined by its stabilizer group S C V n which is an Abelian 
subgroup of V n . To preclude S from containing non-trivial phase factors one usually 
adds a requirement —I ^ S. The set of Pauli operators P G V n that commute with 
all elements of S is called the centralizer of S and is denoted as C{S). If the stabilizer 
group is generated by n — k independent generators, then the centralizer is generated 
by n + k independent generators. Logical operators of a stabilizer code S are elements 
of C(S) which are not in S. One can always choose a set of 2k logical Pauli operators 
Xi, Zi, . . . , X k , Z k e C(«S)\«S obeying the usual Pauli commutation relations. Note that 
C(S) = (i, S, X\, Z\, . . . , Xk, Zk)- Given a stabilizer code S, its codespace is spanned 
by all n-qubit states invariant under the action of S. In this case the codespace is 
isomorphic to a space of k qubits, the so called logical qubits. 

The distance d of a stabilizer code is defined as the minimum weight of a logical 
operator, i.e. 

d= min \P\. (2) 

PeC(S)\s 

A stabilizer code which encodes k logical qubits into n physical qubits and has distance 
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d is denoted as a [[n, k, d}} code. We shall be interested in geometrically-local stabilizer 
codes and similarly local Majorana fermion codes. For such codes, qubits occupy sites 
of a D-dimensional lattice (or, more generally, some graph equipped with a metric) and 
the stabilizer group has a set of geometrically-local generators, S = (Si, . . . , S m ), that 
is, the support of any generator Si has diameter at most r = 0(1) with respect to the 
lattice geometry pj. 

With a geometrically-local stabilizer code S we can associate a Hamiltonian, e.g. 
Hg = — Si where Si is an (over) complete set of stabilizer generators for S. Then the 
ground-space of the Hamiltonian corresponds to the codespace and one may consider 
properties of such physical system at zero or non-zero temperature T. 

Calderbank-Shor-Steane (CSS) codes are a particular subclass of stabilizer codes 
for which the stabilizer group S can be represented as a product of two subgroups, 
S = S(X) ■ S(Z), that contain only X-type and Z-type Pauli operators respectively. 
Any CSS code can be specified by a pair of classical linear codes Cx, Cz Q {0, l} n such 
that 

n 

S(X) = {P = l[X? : (xi,...,x n )eC x }, 
i=i 

and 

n 

S(Z) = {P = l[Z^ : (#]_,..., z n ) G Cz}- 
i=l 

Commutativity between S(X) and S(Z) is equivalent to the mutual orthogonality of 
the codes Cx-.Cz, that is, one must have Y^i=i x i z i = (mod 2) for all x G C x and 
z G Cz- The special subclass of CSS codes obeying Cx = Cz are called weakly self-dual 
CSS codes. For such codes the stabilizer group is invariant under the exchange of X 
and Z operators on every qubit. In that case the code C = Cx = Cz must be a weakly 
self-dual classical code, that is, YH=i x i z i = (mod 2) for all x,z G C . For more 
background on CSS codes, see [22J. 

4. Definition of Majorana fermion codes 

We define a Majorana fermion code as follows. Let Ci, Ci, ■ ■ ■ , c% n be the 2n Majorana 
modes, i.e., operators obeying commutation rules 

The total number of Majorana modes (2n) is always even, because these modes are 
obtained from n original fermionic modes. The single-mode operators ci, . . . ,C2 n , 
together with the phase factor i, generate a group of Majorana operators Maj(2n). 

| Here and below the notation / = 0(g) refers to the limit n — ¥ oo. In particular, r = 0(1) means 
that r is upper bounded by a constant independent of n. 
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Any element of the group Maj(2n) can be represented as r\ ca where A C {1, . . . , 2n} is 
some subset of modes, 

ueA 

and 77 G {±1, ±2} is a phase factor. The subset A is called the support of ca- Here and 
below we use a standard ordering of the product of single- mode operators c u , meaning 
that the indices u increase from the left to the right. We define the weight of a Majorana 
operator as the number of modes in its support, that is, \ca\ — \ A\. A Majorana operator 
is called even (odd) iff its support has even (odd) size. 
For two arbitrary supports A and B we have 

CAC B = (-l) lAm+lAnBl C B C A . (4) 

When either ca or c# is even, the commutation relation only depends on the parity of 
overlap \A D B\. In particular, when regions A and B do not overlap, i.e. \A D B\ = 0, 
then Majorana operators commute, as is trivially the case for Pauli operators on non- 
overlapping supports. However, when ca and cb are both odd and their supports are 
non-overlapping, then ca and Cb anti-commute. 

One can always map the Majorana modes onto Pauli operators on n qubits 
using the Jordan- Wigner transformation T : Maj(2n) — > V n defined as Y(c2i_i) = 
Z\Zi . . . Zi_iXi and T(c 2 j) = Zi . . . Z^Yi, i = 1, . . . , n. Accordingly, the Hilbert space 
Tn describing 2n Majorana modes can be identified with the Hilbert space of n qubits. 
For D > 1 systems, the Jordan- Wigner transformation generically maps local Majorana 
operators onto non-local Pauli operators. This is one of the reasons that local Majorana 
fermion codes defined below may exhibit different properties than local stabilizer codes 
(see e.g. [201 [23] f° r means to make this mapping local by introducing additional modes). 

A Majorana fermion code is determined by its stabilizer group 5 ma j Q Maj(2n) 
which must obey two conditions: 

• <S ma j is an Abelian group not containing —I. 

• All elements of iS ma j have even weight. 

The second condition guarantees that stabilizer operators preserve the parity of the 
number of fermions in the system, and thus any element of 5 ma j is a physically realizable 
operation. Given any pair of stabilizer operators proportional to ca and cb, the 
commutativity caCb = cbCa and the even-weight condition imply that the overlap |y4n5| 
must be even, see Eq. (jlj). 

The set of Majorana operators P G Maj(2n) that commute with all elements 
of <S ma j is called the centralizer of 5 ma j and is denoted as C(5 ma j). Note that the 
centralizer may contain both even and odd Majorana operators. Logical operators 
of a Majorana fermion code <S ma j are elements of C(iS ma j) which are not in iS ma j. If 
<Sma,} is generated by n — k independent generators, then C(<S ma j) is generated by n + k 
independent generators. One can always choose a set of 2k logical Pauli operators 
X 1 , Zi, . . . , X k , Z k G C(<S ma j)\<S ma j obeying the usual Pauli commutation relations. Note 
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that C(<S ma j) = (i, «S ma j, Xi, Zi, . . . , X k , Zjt). In general the set of logical operators 
may contain both even and odd Majorana operators. We give necessary and sufficient 
conditions for a code to contain odd logical operators in Section El see Proposition [TJ 
The codespace of a Majorana fermion code <S ma j is the linear subspace of T n spanned 
by all states invariant under the action of 5 ma j. 

From now on, we shall be interested in geometrically-local Majorana fermion codes. 
For such codes Majorana modes c u occupy sites of some D- dimensional lattice A and the 
stabilizer group <S ma j has a set of geometrically-local generators, iS ma j = (Si, . . . , S m ), 
that is, the support of any generator Si has diameter at most r = 0(1). The codespace 
of 5 ma j coincides with the ground subspace of a Hamiltonian H = — £\ Si that involves 
only geometrically-local interactions among the Majorana modes. 

The simplest example of a local Majorana fermion code is the one describing 
Kitaev's ID model of Eq. (TjQ). The corresponding stabilizer group is 

), (5) 

while the logical operators of the code can be chosen as X = c\ and Z = C2 n - We 
describe a 2D generalization of this code in Section [7] and give some other examples of 
local Majorana fermion codes in Section [HJ 

We can define the distance of a Majorana fermion code similar to the distance of 
stabilizer codes, i.e. as the minimum weight of logical operators, 

d= min |C|. (6) 

C£C(iS ma j)\(S ma j 

According to this definition, a code with distance d is able to detect any error affecting 
less than d Majorana modes, i.e., any operator ca with \A\ < d is a detectable error. 

The notion of a code's distance does not completely capture all aspects of stability 
that Majorana fermion codes can offer since it treats even and odd logical operators 
on the same footing. However, if the system is closed or interacts with a bosonic 
environment, all physically realizable perturbations and error operators must preserve 
fermionic parity and thus must be even. In order to measure the degree of protection 
based on the superselection rules, let us introduce an additional parameter l even defined 
as the minimum diameter of a region that can support an even logical operator, 

^even = mm diam(Supp(C) ) . (7) 

CGC(5 m aj )\<Smaj 

|C|=0 (mod 2) 

As far as zero-temperature stability is concerned, the parameter Z even determines the 
smallest order of perturbation theory at which the ground state degeneracy of the code 
Hamiltonian H = — £V Si can be lifted by a perturbation that involves only even 
geometrically local operators. Hence the parameters d and l even capture two independent 
mechanisms of protection: topological protection by the code distance and protection 
based on the superselection rules respectively. 

Consider as an example the code defined in Eq. (JSJ). Obviously, it has distance 
d — 1. Meanwhile, any even logical operator must include both c±, C2«, and therefore 
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^even = 2n (the lattice size). More generally, using Lemma [2] from Section and the fact 
(see [7]) that any ID stabilizer code has distance 0(1), one can easily show that the 
distance of any ID Majorana fermion code is 0(1). 

It is important to note that the minimum weight logical operator C in Eq. ([6]) 
always has a connected support, that is, one cannot decompose C as C = caCb where 
the separation between A and B is larger than r (the largest diameter of the generators 
of «S maj -). Indeed, in this case both ca and eg would individually commute with 5 ma j. 
Hence ca or cb would be a logical operator which contradicts the minimality of C. On 
the other hand, the minimum weight even logical operator may have highly disconnected 
support as the code in Eq. (E]) demonstrates. 

A simple argument shows that any Majorana fermion code has even logical 
operators. Indeed, if X, Y, Z are logical Pauli operators for some encoded qubit then 
the identity XYZ oc / implies that either all X,Y, Z are even, or two of then are odd 
and the third one is even. We describe properties of even logical operators for general 
.D-dimensional Majorana fermion codes in Section [6j 

As mentioned earlier, the superselection rule prohibiting odd error operators is not 
likely to be completely rigorous. For instance, one might be interested in constructions 
of Majorana fermion codes in which both d and / e ven can be made arbitrarily large by 
increasing lattice dimensions, see Section [7|8l 

Before we continue, let us make a few remarks about general Majorana fermion 
codes based on non-interacting fermions. We would like to make the point that going to 
higher spatial dimensions D > 1 does, in one aspect, not lead to fundamentally different 
behavior as compared to the ID Kitaev's model of unpaired Majorana fermions [15]. 
More specifically, the Bogoliubov transformation allows one to transform any non- 
interacting Majorana fermion Hamiltonian H = i Ylk^i ctkiCkCi into a canonical form 
in which some subset of Majorana modes is unpaired (i.e. these modes do not enter 
into the Hamiltonian). Since each unpaired mode is a linear combination of the original 
Majorana operators c^, the ground subspace of H can be regarded as a quantum code 
with distance d — 1. Hence non-interacting models can only offer protection based on 
superselection rules similar to what the ID model of Ref. [15] achieves. On the other 
hand, unlike in ID, in 2D non-interacting fermion systems with unpaired Majorana 
modes, one can imagine adiabatically changing the Hamiltonian (or 'deforming' the 
quantum code [23]) to move localized unpaired Majorana modes around and enact some 
(but not all) logical gates by braiding. 

5. Code mappings 

In this Section we describe inter-conversions between three classes of codes: (i) qubit 
stabilizer codes, (ii) Majorana fermion codes and (iii) weakly self-dual CSS codes. 

Lemma 1 (Kitaev |25| ) Every [[n,k,d]} qubit stabilizer code S can be mapped onto a 
Majorana fermion code iS ma j on An modes encoding k logical qubits with distance 2d. 
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For completeness, we give the mapping: 

Proof: With every qubit j, we associate four Majorana fermion modes, b*' y,z and 
Cj and hence we have a total of An Majorana fermions. In the n — k independent 
stabilizer generators of S, we replace the local Pauli operators by Xj = iWcj, Yj = ib y Cj, 
Zj = ibjCj. In addition, for each qubit j we add a stabilizer Dj = tfjtfjbjCj to <S ma j (on the 
subspace for which Dj = +1 we have XjYjZj = il). Thus the Majorana fermion code 
<S ma j is generated by 2n — k independent generators, and therefore it encodes k logical 
qubits. Each logical operator of the stabilizer code corresponds to a logical operator of 
the Majorana code. Also, since a logical operator of the Majorana fermion code has 
to commute with each Dj, it must contain an even number of the set {bj,b V j,bj,Cj}, 
and therefore it corresponds to a logical operator of S. Since every Pauli operator 
corresponds to a weight-2 Majorana operator, the distance of the Majorana fermion 
code is twice the distance of the stabilizer code. □. 

Note that by this mapping every operator in the stabilizer iS ma j and logical operator 
in C(5 ma j) will have even weight. In addition, qubit errors get mapped onto even 
Majorana operators. 

Lemma 2 (Doubling) With every Majorana fermion code 5 ma j on 2n Majorana 
modes which encodes k logical qubits and has distance d, we can associate a [[2n, 2k, d]} 
weakly self-dual CSS code. 

Let us first illustrate the idea of the doubling map using the simplest Majorana code 
with 4 Majorana modes and a single generator iS ma j = (C1C2C3C4). Clearly this code has 
k = 1 logical qubit with logical operators X\ = ic\c 2 and Z\ = ic\c 3 . One can easily 
check that iS ma j has distance d = 2. The doubled version of iS ma j is a [[4, 2, 2]] stabilizer 
code with a stabilizer group S = (X1X2X3X4, Z1Z2Z3Z4) obtained by replacing each 
operator c u either with X u or Z u (hence the number of generators is doubled). The 
logical operators of S can be chosen as X\ = X1X2, Z\ = ZiZ 3 , X 2 = XiX 3 , and 
Z2 = Z\Z 2 - 

Proof of Lemma [2j Any operator P G Maj(2n) can be parameterized (up to a phase 
factor) by a binary string x G {0, l} 2n such that multiplication in Maj(2n) corresponds 
to addition of binary strings modulo two. Specifically, one sets x u = 1 if u belongs to the 
support of P and x u = otherwise. Let (ft : Maj(2n) — > {0, l} 2n be the corresponding 
mapping. Consider a classical code 

C = 0(5 maj )c{O,l} 2 ". 

Note that dim (C) = n — k, since 5 ma j has n — k independent generators. Furthermore, 
since iS ma j is an Abelian group containing only even operators, the supports of any 
elements P, Q G <S ma j must have even overlap. Hence C is a weakly self-dual classical 
code, that is, Y^h=i x iUi = (mod 2) for any x,y G C, or, equivalently, 

C C C L . 

Let S = S(X) ■ S(Z) be the weakly self-dual CSS code constructed from C as explained 
in Section [31 By construction, the code S has 2n qubits, n — k independent generators 
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of X-type, and n — k independent generators of Z-type. Hence S encodes 2k qubits. 
Note that each generator of 5 ma j gives rise to a pair of generators in S: the one obtained 
by replacing each single-mode operator Cj with Xj, and the one obtained by replacing 
each single-mode operator Cj with Zj. 

Consider a minimum-weight logical operator in the code S; w.l.o.g. it is either a 
product of X or a product of Zs (but not of both). When we replace each Pauli Xj 
(or Zi) by a Majorana operator Cj, we obtain a logical operator for the Majorana code. 
Vice versa, every logical Majorana operator gives rise to a pair of logical operators for 
the stabilizer code S. Hence the distances of these codes are identical. □. 

Combining the two lemmas, we get the following useful fact. 

Corollary 1 Any [[n, k, d]] stabilizer code can be mapped onto a [[An, 2k, 2d}} weakly 
self-dual CSS code. This mapping preserves geometric locality of a code up to a constant 
factor. 

This result thus shows that in order to derive distance bounds for, say, geometrically- 
local codes, one only needs to prove such bounds for weakly self-dual CSS codes and the 
scaling of rates and relative overhead can be determined by considering only weakly self- 
dual CSS codes. In addition, the code mappings allow one to show that the partition 
function of a Hamiltonian associated with a stabilizer code can be expressed as the 
partition function of a classical Ising (Z 2 ) gauge theory. Since the mapping preserves 
the locality of errors, it is also the physics which is preserved. Hence the properties of Z 2 - 
gauge models, the presence of a phase-transition or not, will be related to the question 
of thermal stability of any stabilizer code [26]. The weakly self-dual character of the 
CSS code, or the fact that the Hamiltonian of the Z 2 -gauge model has only terms with 
even overlap, is crucial. For example, it is well-known that there exists a 3D Ising gauge 
model [27] (in fact, this model is the Z-part (i.e. the subgroup S(Z) of the stabilizer S) 
of the 3D surface code which was shown to be thermally stable against X-errors [28J) 
which has a phase-transition at a non-zero temperature T c . However, this Ising gauge 
model does not have the property that all terms have even overlap. Hence this model is 
not directly pertinent to the thermal stability of 3D stabilizer codes for which one needs 
macroscopic energy barriers against both against X- and Z-error excitations. 

6. Properties of Majorana fermion codes 

In the previous section we have seen that any local Majorana fermion code can be 
mapped to a local weakly self-dual CSS code without changing parameters of the code 
in any significant way, see Lemma [2j Therefore one directly apply any distance upper 
bounds obtained for local stabilizer codes in [7J [29] to obtain analogous upper bounds on 
local Majorana fermion codes. However, one might expect that a Majorana fermion code 
may offer an additional degree of protection resulting from conservation of the fermionic 
parity. Such additional protection may only manifest itself for Majorana fermion codes 
that possess odd logical operators, since for such codes at least some subset of the 
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logical operators is protected by the superselection rules, see Section HI Hence, the first 
question we address in this section is under what conditions a Majorana fermion code 
has at least one odd logical operator. Let us define an operator C a ii measuring the parity 
of the total number of fermions, 

Call = i U ClC 2 ■ ■ ■ C 2n -lC2n- (8) 

Proposition 1 A Majorana fermion code 5 ma j has at least one odd logical operator iff 

Call f. ±<S m aj. 

Proof: Indeed, suppose C a u G ±iS ma j. Then any logical operator P G C (<S ma j ) \£maj 
must have even weight, since it has to commute with C a u. Suppose now that C a n 
±<S ma j. Since the support of C a n has even overlap with the support of any element of 
<5 ma j, we conclude that C a n is a logical operator. Then there must exist another logical 
operator P which anti-commutes with C a u. But this is possible only if P has odd weight. 
□ 

Let us point out that for any Majorana fermion code there is a choice of logical 
Pauli operators such that at most one of them is odd. Indeed, if a code has at least one 
odd logical operator, we can choose logical Pauli operators on the first encoded qubit as 
Zi = C a n an d X\ = P, where P is some odd logical operator, see Proposition [TJ Since 
the logical Pauli operators Xi, Z% on the remaining qubit must commute with Z\, they 
must be even. It will be convenient to introduce a parameter fc odd G {0, 1} such that 
&odd = 1 iff a code has at least one odd logical operator. (It is perhaps important to note 
that fc 0( jd = 1 does not imply that at most one logical qubit has odd logical operators. 
In fact, one can show that there always exists a choice for the logical operators such 
that all logical qubits have logical X and Z operators which are odd weight.) 

For 2D stabilizer codes, it has been proved in [7] that one can always find logical 
operators which are supported on a strip of constant width, i.e., have string-like 
geometry. Lemma [2] immediately shows that the same result holds for 2D Majorana 
fermion codes. In particular, the distance of any 2D Majorana fermion code defined on 
a lattice of size L x L is at most O(L). However, as we mentioned earlier in the paper, 
for Majorana fermion codes with A; dd = 1 we have to focus only on the even logical 
operators while odd logical operators are prohibited by the superselection rules. For 
example, the results of [7J do not rule out the possibility that a 2D Majorana fermion 
code may have logical operators X, Y, Z such that X is even, Y, Z are odd, and X 
has a plane-like geometry, that is, the minimum weight of X is of order n ~ L 2 . Such 
a code might behave similar to the classical 2D ferromagnetic Ising model in terms of 
its thermal stability. Unfortunately, below we prove that 2D Majorana fermion codes 
do not behave like the 2D ferromagnetic Ising model; more precisely we will show the 
following. 

Lemma 3 Let iS ma j = (Si, ■ ■ ■ , S m ) be a local Majorana fermion code defined on a 2D 
lattice (with periodic or open boundary conditions) such that the support of any generator 
S a has diameter at most r — 1 for some constant r. Let A = A! U A 2 U . . . U A t be a 
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partition of the lattice into parallel disjoint strips of width at least r. Then one of the 
following (or both) is true: 

(i) There exists an even logical operator C supported on some strip Ai; 

(ii) There exists a pair of odd logical operators Ci, Cj supported on some pair of strips 
A-i, Aj, ij^j. 

Let us first comment on the implications of the lemma. Consider first the case (i). 
In this case the code has an even logical operator C whose support is confined to a 
rectangular region of size r x L, where L is the lattice size. In other words, at least 
one even logical operator has string-like geometry. Consider now the case (ii). Since 
the logical operators Ci, Cj are odd and have non-overlapping supports, they must anti- 
commute. It implies that Ci and Cj are distinct logical operators, that is, CiCj S ma j. 
But then CiCj is an even logical operator whose support consists of two disjoint string- 
like regions. This result suggests that 2D Majorana fermion codes cannot surpass 2D 
stabilizer codes in terms of their thermal stability by gaining additional protection based 
on the superselection rules. Indeed, as was pointed out by many authors [5J El [30j [7J [9] , 
the existence of string-like logical operators and the lack of string-tension rule out the 
possibility of quantum self-correction at a non-zero temperature. More in particular, a 
logical operator supported on two string-like regions, can be generated by a sequence of 
local even Majorana fermion operators such that for every state obtained in the sequence 
its energy is 0(1) above the ground-state energy. This argument shows that the energy 
barrier between logical states is 0(1), see [7]. 

Note that Lemma [3] can be applied to ID geometry as well by considering a 2D 
lattice of size L x 1. In this situation the strips Ai become intervals of a ID chain of 
length r — 1. Obviously, if a ID code satisfies case (i) of Lemma EJ it does not provide any 
protection at all, since it has an even logical operator of constant weight and constant 
diameter. On the other hand, a ID code satisfying case (ii) of Lemma [3] behaves similar 
to the Kitaev's ID model, see Eq. ([I]). Indeed, such a code has two constant- weight odd 
logical operators Ci, Cj supported on some disjoint intervals Ai, Aj. Clearly, the largest 
possible diameter of the corresponding even logical operator CiCj is of order L. Using 
the notation of Section HJ any ID Majorana code must obey 

d = 0(l) and Z even = 0(L). (9) 

In particular, it shows that Kitaev's ID Majorana chain demonstrates the optimal 
behavior even among the subclass of interacting fermionic models corresponding to 
Majorana fermion codes. 

In Section [7J we construct a 2D Majorana fermion code that demonstrates the 
optimal behavior allowed by Lemma |3j This code has a single even logical operator and 
two odd logical operators with weight of order L located on the opposite boundaries of 
the lattice. 

Lemma [3] can be straightforwardly generalized to any spatial dimension D using 
the partition A = A 1 U A 2 U . . . U A t into disjoint hyper-strips of width at least r, that 
is, rectangles of size s x L x . . . x L, where L is the lattice size and s > r. 
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In the rest of the section we prove Lemma [3J 
Proof: Let us say that a subset of the lattice M C A is cleanable iff for any logical 
operator D G C (<S ma j ) \<S ma j there exists a stabilizer S 1 G <S ma j such that acts trivially 
on M, that is, Supp(Z)S') fl M = 0. Otherwise we shall say that M is uncleanable. We 
shall use the following simple fact. 

Proposition 2 subset M C A zs uncleanable iff' there exists a logical operator 
C G C(5 ma j)\iS ma j supported on M . 

The proposition can be easily proved by applying the Cleaning Lemma from [7] to the 
doubled stabilizer code constructed from 5 ma j as described in Section [5j For the sake of 
completeness, we give a more direct proof of Proposition |2] in Appendix A. 
Now consider two cases: 

(1) There are at least two uncleanable strips Ai, Aj, i ^ j; 

(2) There is at most one uncleanable strip Ai . 

Consider first case (1). Let Ci and Cj be the logical operators of 5 ma j supported on 
Ai and Aj which exist by Proposition [2j If at least one of Ci, Cj is even, we arrive at 
case (i) of the lemma. If both Ci, Cj are odd, we arrive at case (ii). 

Let us now consider case (2). We shall color the strips in black and white in the 
alternating order such that the only uncleanable strip Ai (if any) is black. Then every 
white strip is cleanable. Moreover, since the generators of 5 ma j have diameter smaller 
than the width of a strip, the union of all white strips is also cleanable. As was mentioned 
in Section HI we can always choose at least one even logical operator D G C(5 ma j)\<S ma j. 
Let S G iS ma j be the stabilizer that cleans D from the union of white strips. Then 
C := SD is an even logical operator of S ma j that has support only on black strips. Let 
Ci be the restriction of C onto a black strip Ai. Note that for any i the operator C{ is 
either a stabilizer or a logical operator of <S ma j. If there exists a black strip Ai such that 
Ci is a logical operator with even weight, we arrive at case (i) of the lemma. Otherwise, 
the number of black strips Ai such that Ci is a logical operator with odd weight must 
be non-zero and even (recall that the overall weight of C is even). Hence we can choose 
a pair of black strips Ai, Aj such that Ci and Cj are logical operators with odd weight. 
We arrive at case (ii) of the lemma. □ 

7. Majorana color code 

In this Section we describe a fermionic version of the topological color codes introduced 
by Bombin and Martin-Delgado in [31]. The color codes are weakly self-dual CSS codes 
with geometrically-local generators. One can define a color code on any two-dimensional 
lattice or, more generally, on any surface graph which is 3-valent and admits a 3-coloring 
of its faces. For such graphs two faces share an even number of vertices, as is easily 
checked. Given such a graph, the color code is defined by placing qubits at the vertices of 
the graph. The generators of the stabilizer group are associated with faces of the lattice. 
Specifically, for every face / one defines a pair of generators Sf(X) and Sf(Z) equal to 
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the product of Pauli X's and Z's respectively over all qubits lying on the boundary of /. 
The even overlap condition guarantees that all generators pairwise commute. One can 
show that the logical operators of the color code can be identified with homologically 
non-trivial loops on the lattice, see [31] for details. 

Lemma [2] allows one to identify any color code with a doubled Majorana fermion 
code. A simple example is the 2D color code on a hexagonal lattice with periodic 
boundary conditions (i.e. a torus). Such code encodes 4 logical qubits, see [31]. The 
corresponding Majorana fermion code has a single Majorana mode c u at every site of the 
lattice and a single generator Cf at every hexagon /. The generator Cf is proportional 
to the product of single-mode operators c u over all sites u lying on the boundary of 
/. Note that the pair of generators Sf(X), Sf(Z) can be obtained by applying the 
doubling transformation of Lemma [2] to the generator Cf. Hence the Majorana code 
with stabilizers Cf encodes 2 logical qubits. It is important to note that all 2D color 
codes discussed in Ref. [31] have only even-weight logical operators (with the exception 
of the so-called triangular codes which we discuss in Section E}. The absence of odd- 
weight logical operators in a color code implies the absence of odd logical operators in 
its fermionic version. Hence superselection rules do not play a role in their stability 
properties. 

The formalism developed by Bombin and Martin-Delgado in [31] employs 3-coloring 
of the set of the faces of the lattice to classify the logical operators of the code. As 
we shall see, the global face 3-coloring condition is too restrictive as it leaves many 
interesting color-type codes beyond the scope of the formalism. In particular, the 
Majorana color codes that we describe below are defined on lattices that admit only 
a local 3-coloring, meaning that any topologically trivial region of the lattice admits a 
face 3-coloring but it cannot be extended to the entire lattice. The peculiar feature of 
such codes is that they possess odd logical operators. 

Let E = S 1 x [0, 1] be a two-dimensional cylinder and G C E be a graph embedded 
in E. We shall assume that the graph G induces a cellular decomposition of E, that is, 
the surface E can be decomposed into a set of faces, edges, and vertices that we shall 
denote F, E, and V respectively. The boundary <9E consists of two cycles S 1 x {0} and 
S 1 x {1}. In order to define a Majorana color code, we shall impose four conditions on 
the graph G: 

(Gl) The total number of vertices is even. 

(G2) Each vertex has degree 3 (trivalent graph). 

(G3) The boundary of any face has even length. 

(G4) The boundaries S 1 x {0} and S 1 x {1} have odd length. 

Given a face / G F, let V(f) C V be the set of all vertices that lie on the boundary 
of /. We shall say that a face / G F is adjacent to a vertex u G V iff u G V(f), and 
we say that two faces are adjacent if they share a common edge. The assumption that 
G induces a cellular decomposition of E together with condition (G2) imply that any 
vertex u £ <9E has exactly three adjacent faces and any vertex u G <9E has exactly two 



Majorana Fermion Codes 



16 



adjacent faces. This, together with (G4), implies via Proposition [TJ that we must have 
odd logical operators. We show a non-trivial example of a surface graph G satisfying 
conditions (G1-G4) in Fig. dJ 

Suppose each vertex u G V is occupied by a Majorana mode c u . For any face / G F 
we define a face operator 

c,= n Cu . (io) 

Conditions (G2) and (G3) imply that all face operators have even weight and any pair 
of face operators commute with each other. Using the standard stabilizer formalism one 
can show that there exists a choice of phase factors rjf G {±l,±z}, / G F, such that 
operators rjfCf generate an Abelian subgroup <S ma j(G) C Maj(2n) not containing minus 
identity, i.e. 

SUd(G) = fa/C/, feF). (11) 

is a Majorana fermion code. This code will be referred to as a Majorana color code 
associated with G. One can also regard the codespace of <S ma j (G) as the ground subspace 
of a fermionic local Hamiltonian 

H = -Y.Vf C f ( 12 ) 

(Although the coefficients rjf do not affect any parameters of the code, one needs to have 
an explicit expression for rjf to define the Hamiltonian model Eq. ( fl2l) . Let us mention 
that one can explicitly compute rjf using the geometrical structures known as Kasteleyn 
orientations of a surface [321 [331 E3])- 

We begin by describing the logical operators of the code. Let 7 and 71 be the two 
boundary components of S, that is, 

l0 = V f] (S 1 x {0}) and 71 = V n (S 1 x {1}). (13) 

Define operators 

C Q = J] c u , a = 0,1. (14) 

Note that for any face /, the set V(f) has even overlap with 7 a since one can regard 7 a 
as a boundary of an external face obtained by patching up the hole in S. Therefore C a 
commutes with any face operator. On the other hand, condition (G4) implies that C a 
has odd weight, and thus C a <S ma j(G). We conclude that C a are logical operators of 
the code. In addition, since CqC\ = —CiCq, these are two independent logical operators. 
In other words, the Majorana color code encodes at least one qubit and the logical Pauli 
operators for this qubit can be chosen as 

X = C , Y = C U Z = -iC C 1 . (15) 

The following lemma shows that this is the only logical qubit. 

Lemma 4 The Majorana color code has exactly one logical qubit, i.e., k = 1. 
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Proof: Let us define a classical linear code C C {0, whose codewords describe 
linear dependencies among the face operators. Specifically, a binary string x = {x/}/gf 
is a codeword of C iff 

JJ oc /. (16) 

Using the standard stabilizer formalism one can show that 

k = M - dim (5 maj (G)) = - |F| + dim (C) = dim (C), (17) 

where the last equality follows from the Euler formula |V| + \F\ — = and the 
identity 3|V| = 2\E\. Since we have already shown that k > 1, we know that Eq. ffTB"j) 
has at least one non-trivial solution, that is, 

dim(C)>l. (18) 

Let u be a vertex in <9£, and let / and g be the two adjacent faces. We claim that 
if we fix the values of Xf and x g we uniquely determine a solution to Eq. ffTB"]) . Since we 
see from Eq. ( JIB"]) that x/ ©x 9 = 0, we then proved that there are at most two solutions, 
and hence dimC < 1. 

To prove the claim, let v G V be another vertex of G, and let oj = (iz = 
u,ui, . . . ,u t = v) be any path connecting u and u. Suppose we have already set the 
value of x on some pair of faces adjacent to the vertex it, for some i > 0. Consider 
two cases. Case 1: ^ 9S. Let /ij be the third face adjacent to Ui. From Eq. ( 1TB]) we 
infer that © x 9i © x^ = which uniquely sets Xh v Since two of the faces /j, g>j, /ij are 
adjacent to the edge (w^Wj+i), it sets the value of x on some pair of faces adjacent to 
u i+ i. Case 2: U{ G <9£. In this case U{ has only two adjacent faces fi,gi- If u i+ i dH, 
then both faces fi,gi are adjacent to u i+ i. If G <9£ then only one of the faces fo,gi 
is adjacent to Mj+i, say the face /j. Let = /j and be the two faces adjacent to 
u i+ i. From Eq. ( TIB"]) we infer that © = which sets the value of x on the two 
faces adjacent to Ui + \. This actually even shows that x has to have the same value for 
all faces along a common boundary of E. 

Applying this argument inductively one sets the value of x on some pair of faces 
adjacent to v which sets the value of x on all faces adjacent to v. Since any face is 
adjacent to some vertex, it shows that there is at most one way to extend Xf and x g to 
a solution of Eq. ( 1T61) . □ 

The unique non-trivial solution of Eq. f|T6|) constructed above allows one to define 
subsets of faces F = {/ G F : Xf = 0} and F\ = {/ G F : x/ = 1} such that any 
vertex has exactly two adjacent faces from F\, see Fig. [ljin which the faces from F are 
represented by shaded hexagons. In other words, we have the following corollary. 

Corollary 2 There exists a unique partition of the set of faces F into disjoint subsets 
F and Fx such that each vertex has exactly two adjacent faces from F\ and each vertex 
not lying on the boundary <9£ has exactly one adjacent face from Fq. 
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Another interesting corollary of Lemma H] is that the graph G is not face 3-colorable. 
Recall that a face 3-coloring is a mapping c : F — > {0, 1, —1} such that for any pair of 
adjacent faces f,g one has c(f) ^ c(g). 

Corollary 3 The graph G does not permit face 3-coloring. 

Proof: Indeed, suppose such a 3-coloring exists. Then clearly all faces in F Q must have 
the same color, say, c(f) = for all / G Fq. It implies that all faces adjacent to the 
boundary <9£ must be colored by ±1. However, since the boundary components have 
odd length, such a coloring does not exist. □ 

One can use Corollary [2] to define a family of even- weight logical operators whose 
supports have geometry of a string connecting the two boundary components of E. 
Indeed, define a subset of edges 

E = {e G E : both faces adjacent to e belong to Fx}. (19) 

Any edge e G Eq connects some pair of distinct faces in Fq, or connects some face in 
F Q with one of the two external faces / ex t,o> / e xt,i obtained by patching the holes in X 
(see Fig. [1] where the edges from E are shown in blue). Given any edge (u,v) G Eq 
connecting some pair of faces /, g G Fq, the operator c u c v commutes with any face 
operator Cf, f G F 1: and anticommutes with Cf and C g . Hence we can construct 
logical operators associated with paths of edges in Eq that connect the two external 
faces /ext.o an d / e xt,i- More specifically, consider a graph G^ with a set of vertices 
Fq U / extj0 U / extj i and a set of edges Eq. Let 7 = (ei, . . . , e m ), G E Q , be any path on 
G^ connecting / ex t,o and / ex t,i- Then the operator 

C 7 = H c u c v (20) 

commutes with all face operators Cf, f G F. On the other hand, C 1 anticommutes with 
c and c.\ since it shares exactly one vertex with the boundaries 70 and 71. We conclude 
that C 7 is the logical operator Z ~ cqC\, see Eq. (IT51) . 

The graph G^ defined above allows one to construct a face 3-coloring for any 
topologically trivial region of the lattice. Indeed, as was pointed out above, all faces 
/ G Fq must have the same color, for instance, c(f) = for all / G Fq. Then one needs 
to color the faces / G Fx using the colors c(f) = ±1 such that adjacent faces in Fx have 
different colors. Recall that any pair of adjacent faces in Fx can be identified with some 
edge e E Eq, see Eq. ffl9|) . Hence G admits a face 3-coloring iff the graph dual to G^ 
admits a vertex 2-coloring. Let us denote this dual graph G^\ It has the set of vertices 
Fx and the set of edges Eq. The set of faces of G^ can be identified with Fq. Since 
each face / G Fq has even-length boundary, any homologically trivial cycle in G^ must 
have even length. Hence one can construct a vertex 2-coloring of any subgraph G^ 
that does not contain homologically non-trivial cycles. 

Let us now bound the distance of the Majorana color code focusing on the physically 
relevant case when the generators of the code are geometrically-local. We shall assume 
that S is equipped with a metric such that the boundary components S 1 x {0} and 
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S 1 x {1} have length R, while the distance between them is L. We also assume that 
edges of G have length at most 0(1) and any face consists of 0(1) edges. Below we prove 
that the distance of the code grows linearly with the smallest of the surface dimensions, 
namely, 

d = n(min(R,L)). (21) 
As for the minimum diameter of even logical operators, we prove the bound 

I even = Q\L), (22) 

see Section H] for notations. The regime in which R = 0(1) while L ^> 1 can be 
regarded as protection by the superselection rules only since in this regime the code 
behave similarly to the Kitaev's ID model, see Eq. (J9j). The regime in which both 
dimensions R,L^>1 are of the same order can be regarded as protection by the code 
distance only, since in this regime both even and odd logical operators are equally 
difficult to implement. In the intermediate regimes the code combines both types of 
protection in a way that can be controlled by the choice of R and L. 

Let us now prove the bounds Eqs. (121f22j) . We start from observing that any odd 
logical operator must have weight Q(R). Indeed, any such logical operator P must anti- 
commute with even logical operators 7 constructed above, see Eq. ( 1201) . Obviously, 
one can choose m pairwise disjoint paths 71, . . . ,7 m on the graph connecting the 
two external faces where m = Q(R). Then the support of P must have odd overlap with 
each of the paths ji, i — 1, . . . , m, that is, the weight of P must be at least m. 

Suppose now that P is the minimum-weight operator among the even logical 
operators. Let r = 0(1) be the largest diameter of faces / G F. Consider two cases: 
(i) the support of P can be partitioned into two disjoint components separated by a 
distance greater than r; (ii) such a partition does not exist. In the case (i) we have a 
decomposition P = P\Pi, where Pi, Pi individually commute with any face operator 
and at least one of Pi, Pi is a non-trivial logical operator. Note that Pi, Pi cannot be 
even operators since it would contradict the weight minimality of P. Hence both Pi 
and Pi are non-trivial odd logical operators. However we have already shown that such 
operators must have weight fl(R), that is, we arrive at |P| = |Pi| + |P 2 | = Q(R). Let us 
now consider case (ii). Since P must anti-commute with both c , Ci, its support must 
have odd overlap with both 70 and 71. Condition (ii) then implies that P must have 
weight Q(L/r) = Q(L). In both cases the diameter of the support of P is Q(L) since its 
overlaps with both 7 and 71. 

8. Other constructions of Majorana fermion codes with odd logical 
operators 

In this section we discuss some alternative strategies to construct Majorana fermion 
codes with odd logical operators. Let us begin by considering an unphysical situation 
when the total number of Majorana fermion modes is odd. Then the operator O a n = 
YlueA Cu ^ s an d therefore does not belong to 5 ma j . On the other hand, C a \\ has 
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Figure 1. Example of a surface graph satisfying conditions (G1)-(G4). The lattice has 
periodic boundary conditions along the horizontal axis and open boundary conditions 
along the vertical axis. The subset Fo consists of 10 faces (shaded hexagons). The 
edges of E and E\ are shown using the blue and red color respectively. The boundary 
components 70,71 consist of 5 vertices. 

even overlap with any stabilizer and hence C a n G C(5 ma j)\<S ma j, that is, C a n is an odd 
logical operator. Note however that when the total number of Majorana fermion modes 
is odd, we encode half-integer number of qubits. For example, C a u can be the only 
operator in C (<S ma j ) \<S ma j . Given the mapping between weakly self-dual CSS codes and 
Majorana fermion codes, see Section [5j it is then easy to construct such Majorana 
fermion codes encoding half-integer number of qubits. Indeed, let us take any weakly 
self-dual CSS code [[n, k, d}} where the total number of qubits n is odd. Viewed as a 
Majorana fermion code (see Lemma [2]) it encodes k/2 logical qubits and has n Majorana 
modes. For example, we could take Steane's [[7, 1, 3]] code [22] encoding a single qubit. 
The corresponding Majorana fermion code has three generators, iS ma j = (Si, S 2 , S 3 ), 
where Si = C1C3C5C7, 5*2 = C2C3C6C7, and S3 = c^c^cqCj. The logical X and Z operator 
for the Steane code become a single logical operator (7 a n for the Majorana fermion code 
which encodes half a qubit. Now we can take another copy of this code, or another 
weakly self-dual CSS code with an odd number of qubits and take the product of these 
codes x <S^ aj -. We now have an even number of Majorana fermion modes, hence 
an integer number of encoded qubits. At the same time, odd logical operators of 
and iS^j give rise to odd logical operators of the product code x For two 

copies of the Steane code, the logical X is a weight-3 Majorana fermion operator on the 
7 modes of the first Steane code and the Z is the same operator on the 7 modes of the 
second Steane code. This plug-and-play procedure of adding halves of qubits living on 
separate spatial supports can be enhanced by inserting a piece of passive material which 
encodes no qubits between the two coding regions. The (linear) size of this passive 
region determines the minimum diameter of even logical operators / ev en- In this way 
Kitaev's ID Majorana fermion model can be viewed as a combination of the trivial code 
comprising of a single mode labelled '1', a piece of passive material including modes 2 to 
2L — 1 in which the Majorana fermions are paired, and again a trivial code on mode 2L. 
Such procedure could for example also be applied to another class of 2D color codes, 
namely the triangular codes [31] which encode a single qubit and hence half a qubit 
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when the code is viewed as a Majorana fermion code. 
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Appendix A 

In this Appendix we prove Proposition [2j 

Proof: Given a subset of modes M C A, we shall define two subgroups of 5 ma j. The first 
subgroup denoted as <S ma j (M) contains all elements of 5 ma j whose support is contained 
in M. The second subgroup denoted as contains all operators P G Maj (2n) whose 
support is contained in M that can be extended to some stabilizer. In other words, 
P G <S^j iff Supp(P) C M and PR G <S ma j for some operator R G Maj(2n) such that 
Supp(i?) n M = 0. By definition, one has S maj (M) C S^. 

We shall use the parameterization <fi : Maj(2n) — >■ {0, l} 2n constructed in Section|5l 
see Lemma EJ Consider the linear subspaces (classical codes) 

C = 0(5 maj ), C(M) = 0(5 maj (M)), and C M = <f>(S^). (23) 

By definition, one has the inclusion C(M) C C M . Since the code C is weakly self-dual, 
that is, C C C x , one has XLeA- r «2 /n = ^ ^ or an x e C(M) and y E C. However, since 
a; has support only on M, it translates into ^2 ueM x u y u = 0, that is, we have also an 
inclusion 

C(M) C (C M ) ± . (24) 

By abuse of notations, from now on we shall consider the codes C (M) and C M as linear 
subspaces of {0, l} m , where m = \M\ (note that any vector in C(M) or C M has all zeros 
outside of M). 

There are two possibilities. First, the inclusion Eq. (124]) is an equality, that is, 
C(M) = (C AI ) ± . Taking the orthogonal complement of both sides we get 

C M = C{M) L . (25) 

Let C G C(«S ma j)\<S ma j be any logical operator and x = (f>(0). Decompose x as 
x = Xi Qt © x ex t, where x- in t and x ext have support inside and outside M respectively. 
Since C commutes with any stabilizer supported on M we conclude that x int G C(M) 1 - 
and hence Eq. ( 1251) implies x- int G C M . It means that _1 (xi nt ) can be extended to some 
stabilizer S G <S ma j. Then CS acts trivially on M. Hence M is cleanable. 

The second possibility is that the inclusion Eq. ( |24l) is strict. Then there exists 
some x G 

(C M ) X such that x £ C(M). Let C = Then C has support on M, 

commutes with any element of <S ma j, but does not belong to 5 ma j. Hence C is a logical 
operator supported on M. 



Majorana Fermion Codes 



22 



To summarize, we have shown that if M is uncleanable then the inclusion Eq. (|24p 
must be strict and hence there must exist a logical operator supported on M. 

Let us now prove the converse. Suppose C is a logical operator supported on M. If 
C is odd, then M is uncleanable. Indeed, any stabilizer S G 5 ma j must have even overlap 
with Supp((7) and thus the support of CS contains odd number of modes (and hence at 
least one) from Supp((7) C M. If C is even, then there must exist a logical operator C' 
that anti-commutes with C. Let us show that C'S acts non-trivially on Supp(C') C M 
for any stabilizer S G <S ma j which would imply that M is uncleanable. Indeed, since C is 
even and anti-commutes with C'S, the overlap between Supp(C7) and Supp(C"S') must 
be odd and hence Supp (C"S) contains at least one mode from Supp ((7) C M. Thus in 
both cases M is uncleanable. □ 
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